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ABSTRACT 


The Transonic Turbine Test Rig at the Naval 
Postgraduate School, Monterey, California, uses a super- 
sonic jet exhauster to lower the turbine exit static 
Pressure and thus increase pressure ratios available for 
turbine testing. Exhauster operation results in inter- 
action between two co-flowing streams of air. The primary 
air is supersonic and the secondary air, subsonic. This 
thesis is an experimental study of the mixing process that 
occurs between the two streams in a cylindrical mixing 
Chamber. In addition, a theoretical model of the exhauster 
performance based on polytropic efficiencies and conserva- 
tion of momentum was programmed, and the predicted values 


were compared to test results. 
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Ll. Introduction 

Exhausters are being incorporated into the design of 
many systems because of the simplicity of their design and 
the absence of moving parts. Air augmentation of rocket 
thrust, torpedo drag reduction innovations, use as a jet 
pump in aircraft for fuel transfer or boundary layer con- 
trol are significant applications. 

The Transonic Turbine Tast Rig (TTTR) installed at the 
Propulsion Laboratory of the Naval Postgraduate School uses 
a supersonic exhauster system to increase the range of tur- 
bine pressure ratios available for testing. The test tur- 
bine is enclosed in a hood. A compressor can provide a 
maximum turbine inlet total pressure of 43 psia giving a 


maximum turbine pressure ratio (PE. ' /PS ) of about 3 
nlet exit 


when exhausting to the atmosphere. Since the turbine uses 
only part of the mass flow furnished by the compressor, the 
remainder can be used to drive the exhauster which produces 
a vacuum in the hood and thereby increases the turbine 
pressure ratio. 

Very little information is available on the mixing 
phenomena encountered with supersonic exhausters. 
Abramovich? provides extensive information for the case of 


an incompressible free jet and associated turbulent mixing: 


He includes test work accomplished, and theory developed, 





tapramovich, G. N., The Theory of Turbulent Jets 
(Cambridge, Massachusetts: ‘The MIT Press, Massachusetts 
Institute of Technology, 1963), Part I, "Turbulent Jets 
in Incompressible Fluid.” 
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by Prandtl and Tollmien. The effect of compressibility on 
turbulence is very complicated however, and simplifica- 
tions of the problem must be introduced to develop approx- 
imate solutions. 

This thesis is a study of the supersonic exhauster 
system used with the TTTR. Objectives of the investigation 
were to determine the hood total pressure variation with 
drive nozzle total pressure for a given secondary flow 
rate,and to study the mixing and turbulent interchange that 
occurs between the two co-flowing streams in the mixing 
pipe. Main topics include: traverse probe design, experi- 
mental results of the ejector test series, a computer pre- 
diction program of exhauster performance based on poly-~ 
tropic processes and conservation of momentum, empirical 
determination of polytropic efficiencies, and an attempt to 
determine from test data the losses due to mixing and tur- 
bulent interchange between the two streams of air. 

A summary of the general properties of turbulent jets, 
including terminology used to describe the different mix- 
ing regions that occur with restricted jet operation and 
co-flowing streams, is given in Appendix D. The terms 
ejector and exhauster are used interchangeably throughout 
this thesis. 

The author is most appreciative of the guidance, coun- 
sel and encouragement provided by Professor M. H. Vavra, 


the laboratory assistance given by Mr. J. E. Hammer, and 
the probe fabrication work of Mr. T. B. Dunton. 
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2. Installation. 

The TTTR installation is shown schematically in Fig. l, 
and Fig. 2 is a photograph of the turbine hood and exhauster 
assembly. High pressure air is supplied by an Allis- 
Chalmers VA 312 compressor and enters Tank l. Air to the 
turbine passes from Tank 1 through a calibrated flow mea- 
surement nozzle and enters Tank 2. The amount of air to the 
turbine is controlled by a manual valve; an electrically 
operated exhauster inlet valve controls the drive nozzle 
total pressure (PTN). More detailed descriptions of the 
TTTR and its operating performance are given by Eckert’ and 
by Commons. > 

The exhauster section of the TTTR is shown in greater 
detail in Fig. 3. Important features are the arrangement 
of twelve mixing tube static pressure taps (M2 - M13), the 
Kiel probe to measure drive nozzle total temperature (TTN) 
and total pressure, and the secondary air static pressure 
taps just prior to entry into the mixing pipe. Secondary 
air is defined as the air that enters the mixing pipe from 
the turbine hood. To obtain the mixing pipe static pres- 


Sures, four taps are equally spaced circumferentially at 





2Eckert, R. H., Performance Analysis and Initial Tests 
of a Transonic Turbine Test Rig (Aeronautical Engineer's 
Thesis, U. S. Naval Postgraduate School, Monterey, 
California, May 1966), Section 2, "Installation." 


3commons, P. M., Instrumentation of the Transonic 
Turbine Test Rig to Determine the Performance of Turbine 
Inlet Guide Vanes Through the Application of the Momentum 
and Moment of Momentum Equations (Master's Thesis, U. S. 
Naval Postgraduate School, Monterey, California, 1967), 
Section 2, "Installation." 
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each location and commoned to obtain an average reading. 
Secondary air static pressures are measured by three taps 
at each location, equally spaced circumferentially and 
commoned. 

One section of the mixing pipe was modified by cutting 
a 3 inch by 1/4 inch slit at 3 locations to provide ingress 
and egress of a traverse probe. The slits were designed to 
Place the traverse probe total pressure opening at the same 
distance from the drive nozzle exit as the mixing pipe wall 
static pressure; that is, in the plane of the wall static 
pressure taps for a given axial location. Seals were also 
fabricated to prevent leakages, with or without the tra- 
verse probe in position. 

Secondary air total temperature (TTD) is measured with 
a modified Kiel probe. Hood total pressure (PTD) is mea- 
sured by eight static pressure ports that are commoned to 
obtain an average value. The secondary air is essentially 
that which passes from Tank 2 to the turbine. It is aug- 
mented by a small amount that enters the shaft labyrinth 
seal of the hood. Dimensions of the shaft labyrinth shown 


4 due to increas- 


in Fig. 4 differ from those given by Eckert 
ing shaft clearance from 0.005 inch to 0.020 inch. A 
thermocouple is placed at the labyrinth entrance to measure 
the total temperature (TTC) of the leakage air. 


Drive nozzle I, the exhauster drive nozzle designed for 


the TTTR, is shown in Fig. 5 with detailed measurements. 


4Eckert, op. Git. , Pay, p. 25. 
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Drive nozzle II, a smaller nozzle with the same area ratio 
(Aguee0/ME ay) was designed and built to investigate the 
effects of different primary and secondary flow areas, and 
smaller primary flow rate,on the mixing process and on hood 
total pressure. Details of Drive nozzle II are shown in 
Fig. 6. Design calculations were made to insure a uniformly 
decreasing secondary flow area for the new nozzle. 

All temperatures are measured using standard Iron - 
Constantan or Cromel - Alumel thermocouples with an ice bath 
as reference. Differential pressures across the flow mea- 
Suring nozzle are measured in inches of water on a mano- 
meter board with 0.1 inch gradations. Mixing tube static 
pressures are measured against atmospheric pressure with a 
Texas Instruments Fused Quartz Pressure Gage. All other 
pressures are measured in inches of mercury against atmos- 
Pheric pressure as reference using 96 inch long manometer 
tubes with 0.1 inch gradations. A check of static pressure 
readings between the Texas Instrument gage and the mercury 
manometers produced values within 0.5 per cent of each 


other. 


3. Probe Design. 





In order to study the mixing of two co-flowing streams, 
velocity traverses were desired at a number of axial loca- 


tions. Naviaux> attempted to obtain this information but 


"Naviaux, J. C., Transonic Turbine Test Rig Exhauster 
= Tests and Tests of a Reaction Turbine (Master's 
esis, U. S. Naval Oca. aed School, Monterey, 
California, December 1966), Section 4, "Test Results." 


vy 


experienced difficulty in obtaining reliable data. During 
his tests the probe had excessive vibration, the static 
Ppresgure probe used had the pressure ports only 2 diameters 
away from the probe tip, and a thermocouple was not incor- 
porated, hence, the differences in stagnation temperatures 
between the two streams could not be measured. 

For the ejector tests, the four probes of Fig. 7 were 
designed and built to obtain traverse information; two were 
pressure probes with a thermocouple and two were hot wire 
probes. Probe I was the only successful instrument and all 
traverse information was obtained with it. 

Probe I is shown with its significant dimensions in 
Fig. 8. Criteria presented by Gettelman and Krause® con- 
cerning static port location for maximum accuracy were 
strictly adhered to. Best results are obtained when the 
static ports are 4 - 6 sensing element diameters away from 
the tip,and greater than 8 support diameters away from the 
support. These minimum criteria produce static pressure 
readings within 5.0 per cent of the velocity head at 
Mach 0.9 and 3.5 per cent at Mach 0.5. Mr. E. Keener of 
NASA Ames Laboratories assured the author that the error 


decreases in supersonic conditions. 





Scettelman, Cc. Cc. and Krause, L. N., Considerations 
Entering into the Selection of Probes for Pressure 


Measurement in Jet Engines (ISA Proceedings Vol. 7, 
Aeronautic Instrumentation - Session II), pp. 134-137. 
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The blockage of probe I is 4.3 per cent of the mixing 
Pipe area. mata" show that a blockage effect based on ratio 
Of support diameter to duct or jet diameter can induce an 
error as large as 8.0 per cent for conditions and dimensions 
used in the ejector tests. The error, however, causes a 
measurement that is too large in a circular duct anda 


measurement that is too small in a jet. This difference in 
direction of the measurement error probably leads to a 


partial cancellation of the blockage error at mixing stations 
near the drive nozzle for the exhauster system. The exact 
effect and variations with axial position were not evaluated. 


The thermocouple of probe I was placed in the support 


0.5 inch above the pressure element. This, position was 


chosen to avoid a second extended sensing arm which could 
induce errors due to probe interactions. The thermocouple 
is exposed to the stream by a 0.160 inch diameter hole. A 


0.040 inch diameter hole in the after side provides air 


circulation. Total temperature measurements at mixing pipe 
station Ml were less than 2.0 per cent different than those 
obtained for the primary and secondary flow before entry to 
the mixing pipe. 

The probe vibration problem was eliminated by instal- 
lation of two 4 inch x 6 inch box supports contoured to fit 


the outside surface of the mixing pipe. Each box was 


4 
Krause, L. N. and Gettelman, C. C., Effect of Inter- 

action Among Probes, Supports, Ducts, Walls and Jet 

Boundaries on Pressure Measurements in Ducts and Jets (ISA 


Proceedings Vol. 7, Aeronautic Instrumentation - Session II), 
Fig. 4, p. 140. 


fitted with a 3 inch teflon bushing to provide probe support 
over a long distance. This eliminated the pivot point about 
which a probe would tend to oscillate. MThe assembly is 
shown in detail in Fig. 9 and pictured in Fig. 10. The 
design proved very successful. 

The mixing process is primarily a function of axial 
location in the mixing pipe. Probe I has a 3 inch axial 
distance between the position where the total pressure is 
taken and the thermocouple. Because significant gradients 
in total temperature exist, particularly near the drive 
nozZle, probe II was designed to obtain the temperature 
nearer to the total pressure opening. An exposed thermo- 
couple was placed in a position 1/2 inch aft of the static 
pressure ports. 

Testing of probe II produced expected results in the 
subsonic secondary flow. In the supersonic region of the 
primary flow, the test data were irregular both for static 
pressures and temperatures. However, lower temperatures 
were recorded whereas in actuality the primary flow had a 
higher total temperature than the secondary flow. The lower 
temperature reading occurred because an exposed thermocouple 
in supersonic flow measures an adiabatic wall temperature. 


Shapiro® introduces the energy equation to describe the 


8 
Shapiro, A. H., The Dynamics and Thermodynamics of 


Compressible Fluid Flow Volume II (New York: The Ronald 
Press Company, 1954), pp. 1034-1036. 
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relationship between the adiabatic wall or recovery tempera- 
ture Tay, the static temperature T,, and velocity V 
that must exist in the boundary layer over such an exposed 


thermocouple. 


2. 
Te = te +t RE, 


R is the recovery factor that has been shown in many experi- 


2 to be equal to the square root of the Prandtl number 


ments 
or 0.85 in laminar flow. Probe II could be used in super- 
sonic flows only if the velocity were Known and actual tem- 
perature, the quantity to be determined from measurement. 
This temperature was needed in the ejector tests to deter- 
mine velocity so probe II was not applicable. 

An attempt was made to verify the shape of the pecul- 
iar velocity profiles obtained at axial stations Ml - M3 
by using a hot wire anemometer. The results of such a test 
would be qualitative because of the temperature difference 
between the primary and secondary streams. However, the 
difference in total temperature never exceeded 80°F. With 
the hot wire temperature at 1200 OF, results within 10 per 
cent accuracy could still be obtained. This would be suffi- 


ciently accurate to confirm the velocity profile obtained 


with probe I. 


Probe III is shown in Fig. 7. It is a hot wire mount- 


ed on vertical arms attached to 1/4 inch tubing. On each 


9 
Ibid., pp. 1070 - 1072. 


call 


of six tests the wire maintained its integrity in the secon- 
dary flow but broke as soon as the primary jet was en- 
countered. It is possible that the hot wire broke because 
Of vibrations of the thin vertical arms to which the wires 
were attached. 

Probe IV, shown in Fig. 7, was designed such that the 
arms to which the wire was attached were stronger and ori- 
ented parallel to the flow. Both changes were made to 
dampen arm vibration. The probe was placed in the center 
stream before the exhauster was operated. The air flows 
through the turbine and drive nozzle were slowly increased, 
but the wire broke at velocities far below those experi- 


enced during normal exhauster operating conditions. 


4. Description of Test Runs. 

Tests were made to determine the hood total pressure 
Variation as a function of drive nozzle total pressure and 
secondary flow rate,to establish the exhauster performance, 
to obtain the variation of mixing pipe static pressure as a 
function of drive nozzle total pressure and secondary flow 
rate, and to study the interaction of the two co-flowing 
streams by obtaining velocity profiles at different axial 
positions of the mixing pipe. Testing was accomplished 
with two different drive nozzles in order to compare the 
effect of variation of the primary flow rate and secondary 
flow annulus area on the problem. 


A total of 34 runs was carried out between February and 


September 1967. Exhauster performance and mixing pipe 
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static pressure rise were obtained for drive nozzle total 
pressure ratios (PIN/P tm) Of 2e65ye 2245, -2eb6, 1a90m amd 
1.63 and for secondary flow rates of 4 and 2 lbm/sec. 

The TTTR was set in the following manner. Equation 1 
given in Sec. 5,entitled "Data Reduction”, is the equation 
used to determine actual flow rate through the flow measur- 
ing nozzle. Using average values for co. # Kpeand Yy 
and with D> known for the TITR installation, an expression 
for the differential pressure h, across the pressure taps 


Was obtained. 


alialsae 
h = Lae VE 


Ww 
w = desired secondary flow rate 

Py = static pressure at upstream tap 
Ty = fluid temperature at upstream tap 


For each test condition, the drive nozzle total pressure was 
set by hand. After allowing time for stabilization of con- 
ditions, readings of the calibrated flow measuring nozzle 
total temperature Ty and pressure P,; were taken,and an 
hy, calculated that corresponds to the desired flow rate. 
By adjusting the turbine inlet valve, this value of h, 
was set. Flow rate regulation to within 0.05 lbm/sec of 
the desired amount was obtained. 

There was no attempt made to precisely regulate either 


primary or secondary air temperature due to the large number 


of variables and valve settings involved. Turbine testing 
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being carried out concurrent with the ejector tests created 
different turbine efficiencies and operating conditions 
which made temperature variations difficult to regulate. 

The analysis of early test results showed that traverse 
probe position significantly affected hood total pressure 
and secondary flow polytropic efficiencies. The probe was 
removed during final tests to obtain exhauster performance 
and mixing pipe static pressure rise. 

Velocity traverses were initially made at mixing sta- 
tions Ml, M2, M3, M4, M6, M7, and M8 of Fig. 3 for a pres- 
sure ratio PTN/P.,,, Of 2.18,and secondary flow rates w, 
of 4, 3, and 2 lbm/sec. Data points were taken in 1/2 inch 
radial increments at stations M6 —- M8. Due to large gradi- 
ents in total temperature and pressure, the increment was 
changed to 1/4 inch for stations M2 - M4 and to 1/10 inch 
at Ml. Additional traverses were obtained at Ml and M4 for 
w, of 4 Ilbm/sec and PIN/P,,, of 2.65 and 1.63. 

It was not possible to change the traverse probe posi- 
£ion without shutting down the exhauster 5 Sea After 
restart and stabilization of the system, the new operating 
conditions could vary in total temperature up to 3 degrees 
from the previous run. The consequence was considered 
secondary compared to the effect probe position had on ex- 
hauster operation. 

Testing with probe II in run 17 produced total pres- 
sures different from those obtained with probe I for the 


same operating conditions. Investigation of the equipment 
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showed a leak in the traverse probe total pressure tubing. 
The effect of the leak became noticeable with probe II 
because Of its smaller total pressure opening. After re- 
pairs, a retest of data points indicated the leak had 
existed during all traverse tests. 

The planned test program was revised and retests ini- 
tiated. Final tests included obtaining the hood total 
pressure and mixing pipe static pressure rise versus drive 
nozzle total pressure, for secondary flow rates of 4 and 2 
lbm/sec. Velocity traverses were taken at stations Ml - M8 


(less M5) for w,. of 4 and PTN/P of 2.18. The tests 


t 


were carried out for both drive nozzles. 


atm 


Several attempts to verify the velocity profile at 
station Ml with a hot wire anemometer proved futile. The 
wire could not withstand the impact of the primary stream 


and broke immediately upon insertion in every test. 
5. Data Reduction. 
Secondary Flow 


Secondary flow rates were measured by a flow nozzle 


between tanks 1 and 2. ‘The basic formula for flow rate in 


lbm/sec was given by Eckert 19 to be 
e Zz - 
Ww = OG Dacca a ity - 


1Onckert, op. cit., Eq. 1457 p. 98. 
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D5, = nozzle throat diameter (in.) 


OC = coefficient of thermal expansion based on 


Ds (dimensionless) 


K = discharge coefficient (dimensionless) 

Y, = expansion coefficient (dimensionless) 

P, = static pressure at upstream tap (in. Hg.) 
T, = fluid temperature at upstream tap (°R) 


h, = differential pressure across taps (in. H50) 


Values for C€ were determined as outlined by Eckert. 11 


The expansion factor Yj was changed to that given by the 


ASME Power Test Codes.12 Discharge coefficient K was 


13 to be a function of Reynolds number. 


found by Naviaux 
Naviaux plotted nozzle flow coefficient versus Reynolds 
number and by curve fitting obtained an eighth order poly- 


nomial for K, with Reynolds number the independent varia-. 
ble. Computation of the polynomial showed that even the 
highest order terms produce large corrections for K. For 


Re Of 11.31 xX 10° the seventh order correction was +1168.857 


tlipid., pp. 98-99. 


i 
Flow Measurement, Supplement to ASME Power Test 
Codes, New York: ASME, 1959, Chap. 4, Part 5, p. 74. 


13Naviaux, Ge. Cilt., Pp. So. 


26 


and the eighth order term -96.291. Because values of K 
were determined by Naviaux to be in the small range from 
1.00 to 1.05, the polynomial with such large corrections 
was considered inappropriate, and K was varied in steps 
of 0.01 as a function of R, using Naviaux's graphical 
14 


results as a reference. 


The method used for estimating the leak rate through 


the shaft labyrinth was developed by Eg1it° and revised by 
Eckert l® to apply to the TTTR installation. 
. AR, \9 
West = a P _ i, (2) 
a 
S Te, R 
and 


1- (%,,) 
a 


n+ In("e/-p) 





(3) 


CC = discharge coefficient for a single 
throttling (dimensionless) 

oc = 0.76 from Egli for present installation 

A = cross sectional area of seal passage (in. *) 


l4tnid., Fig. 19, p. 68. 


Seg li, A. "The Leakage of Steam Through Labyrinth 
Seals," Transactions of ASME, V. 57, 1935, pp. 115 - 122. 


‘nckert, op. cit., pp. 103 - 105. 
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PL = total pressure at seal entrance (psia) 
Pp = static pressure at seal exit (psia) 

Te = total temperature at seal entrance (°R) 
n = number of throttlings 


The calculation of secondary flow rate was accomplished 


with computer program SECFLO which is shown in Appendix C. 


Velocity 


The velocity profile data reduction was accomplished 
by the computer program VELOCITY. Total temperature, total 
pressure and static pressure were measured. The pressure 
ratio was tested for supersonic or subsonic conditions. 

17 


For supersonic flow, the Rayleigh Pitot tube formula was 


used. 


(4) 





For subsonic conditions, the standard compressible flow 


relation was used. 


Boz (j+ stm?) (s 





17 snapiro, A. H., The Dynamics and Thermodynamics of 


Compressible Fluid Flow Volume I (New York: The Ronald 
Press Company, 1954), p. 118. 
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In either case, Mach number was solved by iteration and 


converted to velocity with the measured total temperature, 


using: 


Tes jy X=! M* => ilee= 


aim R. | + Ley Mm @ (6) 
2 
Cel, = Cel +5 => V= \2qIee(E-T) oa, 


The computer program used for the reduction is given in 


Appendix C. 


6. Exhauster Performance Prediction. 

An analysis of the exhauster performance to predict 
the hood total pressure, for specified drive nozzle total 
temperature and total pressure, and secondary air mass flow 
rate and total temperature, has been carried out by Vavral8 
in the TTTR preliminary design calculations, and used by 
<— However, both analyses assume that the drive 
nozzle is operating at design conditions, and the static 
pressure at the entrance to the mixing section is constant 


in radial direction. Considering the many variables in the 


eR M. H., Theoretical Evaluation of Exhauster 
System of Transonic Turbine Test Rig (TN65T, U.S. Naval 
Postgraduate School, Monterey, California, 1965), pp. 16 - 


21. 


Way iaux, op. Cit., pp. “22 - Za 
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physical set-up of the TTTR, such as turbine efficiency, 
RPM, secondary air mass flow rate and total temperature, 
and the range of drive nozzle total pressures available, it 
seems unlikely that the nozzle will operate always at de- 
Sign conditions. Test data plotted in Fig. 29 indicate a 
static pressure gradient at station Ml, located just 9/16 
inch from the nozzle exit. Therefore, the possibility of 

a discontinuity in static pressure between the nozzle exit 
(area Al) and the secondary flow annulus (area A1T) at the 
start of the mixing section is contained in the present 
development. This hypothesis is supported by Jaumotte*? 
who indicated that the classic hypothesis that assumes 
uniformity of the static pressure in the initial section of 
mixing is wrong. In the cited reference it is stated also 
that the laws governing the pressures in the two fluids and 
the factors which influence them are not yet well under- 
stood. 

The complete revised analysis of exhauster performance 
is included as Appendix A. Appendix B gives the computer 
program used to calculate prediction maps of exhauster per- 
formance over its operating range. The pertinent steps of 
the analysis are outlined below with the notations corres- 


ponding to those of Fig. 3. 


20mngle, M. O., Some Problems in the Design and 
Operation of Jet Ejectors, Proceedings of the Institution 


of Mechanical Engineers Steam Plant Group (London, England, 
Ne. 33, 968), p.350. 


Polytropic process relations as well as non-dimensional 
Zi 
flow functions are used. These are taken from Vavra. For 


an expansion process 


(8) 














where 
nN ‘s 
¥- 2, (¥-1) 
n = polytropic exponent 
; - at 
Np = polytropic efficiency aT, 
6 = ratio of specific heats 
w = flow rate (lbm/sec) 
aS = inlet stagnation temperature (oR 
O 
Ag = flow area at station of interest (in.°) 
P, = inlet stagnation pressure (psia) 
O 
P, = static pressure at station of interest 
(psia) 
g = universal gravitation constant 


(lbm-ft/lb-sec?) 
Subscript 1 refers to inlet conditions 


Subscript 2 refers to station of interest 


elvavra, M. H., Problems of Fluid Mechanics in Radial 
Turbomachines (Rhode-Saint-Genese, Belgium: Von Karman 
Institute for Fluid Dynamics, 1965), VKI Course Note 55a, 
pp. 22 - 35. 


31 


For a compression process 








eres 
h\w 


iy (9) 





ay! le, / th 


alee) 


P, = static pressure at 1 (psia) 
ep 
Ls e(i-%,) 


For a choked nozzle 





v IT, 
cae * o, 


0 


nl = 


ia 

Ni-] 
$ J [3 
~ choked N+! 


Initially, polytropic efficiencies were assumed to be 


(10) 





constant over the operating range of the exhauster. The 
determination or estimation of the values used is described 


in Sec. 8 entitled "Determination of Polytropic Exponents.” 


Drive Nozzle 


The flow functions from the nozzle entrance to the 


throat ( ¢, ), and from the nozzle entrance to the exit 


( $, ), are 


e = i ae ya (iy) 


Wr LT R sia 
J 


: Ay FN 


e, 


a2 


The relationship 


oA, = oA, 


(13) 


exists,and 


A 
$= 2n B, ~- 


Expressing the exit flow function in terms of a pres- 


sure ratio 












$s 


2x 


Palle a 


With P .., given by Eq. 


nt 
| Y) (15) 
Pen 


14 and PTN an initial condi- 


tion, P1N may be solved by iteration. With PIN, the 


quantities T1N, VIN, and we may then be obtained. 


Secondary Flow 


For the expansion process from the hood to the mixing 


Pipe annulus area, there is 


/-—o (z ent fe Wes 
man = > = s — ( 3 ek (16) 
(=) Tro 


With secondary flow rate and hood total temperature 


known, a value of P1T is assumed and PTD calculated by 


iteration and matching. For the assumed P1T and correspond- 


ing PTD, the values of T1T and V1T can be calculated 
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Mixing Pipe 
With the values calculated for the drive nozzle and for 
assumed secondary air flow conditions, the momentum equation 
was used for the region within the mixing pipe to obtain 


conditions at the diffuser entrance. The basic equation is 


( Wet), - ¥eV,, ~My, = (PIT)AIT + (PIN)A] ~(P2)A2-F, ayy 
Assumptions include uniform velocity profiles, uniform sta- 
tic pressure distribution over each area, and adiabatic 
flow. The term Fr, represents the viscous shear force 
imposed by the wall on the moving fluid. Estimation of the 
magnitude of this force is outlined by Naviaux. In the 
calculation, the resultant pressure drop produced by Fe 
is determined by a single iteration of Eq. 17. 

In Eq. 17, Py and V2 are unknown. With continuity, 
the pressure PP) can be expressed with the known quanti- 
ties and V5, which leaves one equation to be solved with 
the unknown quantity V5. With V5, all conditions at the 


diffuser entrance can be obtained. 


Diffuser 


For the compression in the subsonic diffuser 


ae fs E e \* — ail 
ae a ; 
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The pressure P3 is ebtadinedeby eitemation; and) if) it dees 
not match atmospheric conditions, a new value of PIT for the 


secondary air must be assumed and the calculation repeated. 


7. Test Results and Discussion. 


Hood Total Pressure 

Tests of the exhauster system were conducted to measure 
the hood total pressure obtained aver the operating range of 
the drive nozzle. Runs were made at a secondary flow rate 
of 4 and 2 lbm/sec for each drive nozzle. 

Two attempts to obtain the information were discarded. 
Run 11 was made with the traverse probe at station M2 of 
the mixing pipe. Comparison with other test data showed 
that hood total pressure was increased by as much as 1 inch 
of mercury with the traverse probe installed, particularly 
when the probe was close to the drive nozzle. The increase 
in hood pressure is the result of losses due to shock waves 
that attach to the probe, the drag of the probe in the 
stream, and the wake formed behind the probe. Run 29 with 
drive nozzle II was not used due to the discovery that a 
large seal in the hood had not been properly installed. 
Run 18 and 22 for drive nozzle I and run 32 for drive nozzle 
II were successful. The results are plotted in Fig. ll. 

The graph shows that for the same secondary flow rate, 
drive nozzle I creates a lower hood total pressure than 
nozzle II. This occurs primarily because nozzle I has a 


2 


throat area of 10.42 in.“ and nozzle II a throat area of 
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8.29 in.2 so that with choked conditions, about 1.51 lbm/sec 
greater flow rate of primary air exists with nozzle I. 
Hence, its driving momentum flow is larger. 

For drive nozzle I, a lower hood pressure resulted for 
a smaller secondary flow rate at a given drive nozzle pres- 
sure ratio. In a relative sense, the effect is the same as 
described above. For the same drive nozzle pressure, which 
establishes the mass flow rate and momentum flow of the 
primary nozzle, it is possible to produce lower pressures 


in the hood for smaller secondary flow rates. 


Mixing Pipe Static Pressure Rise 
Mixing pipe static pressures were recorded for each 


pressure ratio PTN/P and flow rate. Graphs showing 


atm 
the static pressure variation as a function of distance 
from the drive nozzle exit are given in Figs. 12, 13, and 
14. In general, there occurs a rapid increase in the sta- 
tic pressure along the initial region of the mixing pipe. 
In the main mixing zone where the primary and secondary 
core regions have disappeared, there is a slight but steady 
decrease in static pressure in each case, except for 


PTN/P = 2.65 in Fig. 12. The exception could be due to 


atm 
the fact that incomplete mixing has occurred in the mixing 
pipe for the extreme condition tested. This point is the 


greatest drive nozzle pressure ratio used and lowest hood 


total pressure produced during testing. As expected, the 
static pressure rises in the diffuser to the atmospheric 


pressure at the exit. 
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No attempt was made to predict analytically the mixing 
pipe static pressure as a function of axial location. Due 
to the extremely complex nature of the flow which includes 
mixing, turbulence, and shocks in this region, the momentum 
equation used in the exhauster prediction analysis is 
applied only for the whole mixing pipe to obtain average 
conditions ahead of the diffuser. The rapid rise of static 
pressure in the initial region is due to the slowing down 
of the supersonic stream as it mixes with the subsonic 
secondary air. Also present is a series of expansions and 
shock waves to equalize the static pressure between the two 
flows. A possible extension of this thesis would be to 
devise a mixing chamber for Schlieren photographic study 
of the mixing phenomena in compressible flow. 

The slight but steady decrease in static pressure in 
the main mixing zone indicates the loss of available energy 
due to losses caused by the shear forces in the flow and 
the continued turbulent interchange between particles. 

Figures 13 and 14 show that at L/D of about 4, the 
wall static pressure is higher for a drive nozzle total 
pressure ratio of 2.45 than for 2.18. This condition re- 
peated itself on 4 different occasions, but each time the 
manometer board was very unstable with pressure variations 
at stations Ml — M3 up to plus or minus 3 in. Hg. observed 
out of 6 - 8 in. Hg. vacuum. The unsteady condition was 
not reflected in the hood total pressure. It appeared to 


be the result of a shock moving fore and aft in the initial 


37 


region of the mixing pipe. A very definite sound variation 
was heard while operating at this condition which would 
tend to support the hypothesis of an unstable shock move- 


ment. 


Velocity Profiles 

Velocity profiles obtained with probe I for drive 
nozzle I are shown in Figs. 15 through 19, and those for 
drive nozzle II in Fig. 20. There occur large gradients 
between the two co-flowing streams at station Ml and con- 
siderable non-uniformities of the velocity profiles. Both 
of these effects decrease with increase in axial distance 
from the drive nozzle. 

The non-uniformity of velocity profile had been ex- 
pected because of Naviaux's results? but not to the degree 
observed in the ejector tests. However, Naviaux did not 
traverse closer than 4 inches to the drive nozzle, whereas 
at station Ml the traverse total pressure opening is only 
9/16 inch from the drive nozzle exit plane. 

The unusual profiles are probably caused by the shock 
structure that actually exists in the flow. A full discus- 
sion with sketches of the possible pattern is included in 
Sec. 10. Essentially, the presence of internal oblique and 
normal shocks causes sharp gradients in pressure and veloc- 
ity in the flow with accompanied differences in flow direc- 


tions. The large peak in static pressure shown in Fig. 29 


22 naviaux, op. cit., Fig. 3) pewoe. 
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for station Ml is probably due to a large radial velocity 
component in the mixing zone at this radial position. 


23 for a free 


Abramovich cites experiments by Zhestkov 
jet operating at off-design conditions. The results are 
similar to the ejector tests with regard to the decrease 
in velocity in the center of the primary jet,compared to 
that in the jet’s boundaries. By comparing the velocity 
profiles at stations Ml, M2, and M3 it is apparent that the 
secondary air increases in velocity and the secondary core 
region decreases in size. The primary stream velocity 
decreases and the primary core area decreases as the mixing 
zone increases in size. 

Interesting results were obtained by non- 
dimensionalizing the velocity profiles obtained in the 
Main mixing region of the mixing pipe,in the manner given 
by Abramovich?“ for compressible flows. The main mixing 
region is the part of the mixing pipe where the "core" 
flows of primary and secondary air have jueoemeceade | 
Appendix D describes the regions more fully. Abramovich 
uses the ordinate V/V, where V is the velocity at a given 


location, V is the axial velocity. For the absissa, 


ax 
sy 
§ = FT 


23 .hramévich, “op/ eit., Bp. 335 - 338. 


2rd, pp. 0275 eeaaee 


oo 


where y and ya represent the coordinate of the point 
which corresponds to an arbitrary value of the dimension- 
less velocity,and the coordinate where the velocity is 1/2 
the axial velocity,respectfully. Figures 21 and 22 graph 
the result. In both cases a single profile is obtained. 
Station M4 is not inside the main mixing region with 
drive nozzle I,but it is with drive nozzle II. Figure 23 
shows that the non-dimensional profile at station M4 for 
nozzle I deviates from the single profile obtained by non- 
dimensionalizing the velocity profiles of stations Mo —- M8. 
From these results, it appears that a sufficient condition 
to determine the start of the main mixing region in a 
restricted jet is to determine where the velocity profile 
may be non-dimensionalized to a single profile. The 
following paragraphs discuss the mixing regions in greater 


detail. 


Mixing Zone 
The shape of the mixing zone determined from the velo- 
city profiles is shown in Figs. 24 and 25 for nozzle I and 
nozzle II, respectively. This is at best an approximation 
because of the limited number of traverse positions in the 
initial mixing region. The mixing zone appears to touch 
the mixing pipe wall at an L/D of 3.78 for nozzle I,and 


4.71 for nozzle II. D is the nozzle exit diameter. 
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Victorin*> Carried out a series of experiments using 
water as the driven and driving fluid and determined that 
the velocity distribution in the mixing chamber behaved in 
the same fashion as a circular wake with the width propor- 
tional to 3 where x is the axial position. The 
assumption is made that the primary fluid emanates at a 
point source at x = 0 located at a position upstream of 
the nozzle so that the velocity distribution would have 
spread to match the nozzle exit dimensions. Calculations 
of the theoretical distance xl/3 | to predict where the 
mixing zone meets the mixing pipe wall, gave an L/D of 8.27 
for nozzle I,and 9.71 for nozzle II. These conditions are 
shown in Fig. 26. 

It appears from the traverse information that the out- 
ward spread of the mixing zone occurs at a faster rate than 
the incompressible approximation. From the test data ob- 
tained, a closer approximation for the spread of a jet in 
compressible flow would be y = res The relation predicts 
that mixing zone touches the pipe wall at an L/D of 3.87 
for nozzle I,and 4.86 for nozzle ITI. 

From the traverse data from guns 4 - 17 it can be con- 
cluded that the shape of the mixing zone is independent of 


secondary air mass flow rate. This conclusion is valid 


even though a leak existed in the traverse probe total 


25victorin, K., Investigation of Turbulent Mixing 
Processes, NACA T 10 96, 1946. 
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pressure line during these runs. The leak did not affect 
the radial position of the steep total pressure gradient 
observed in test data. The radial position where the steep 
gradient in total pressure occurred was invariant as the 


secondary flow rate was set at 4, 3, and 2 lbm/sec. 


8. Determination of Polytropic Exponents. 


General 
The theoretical model which has been assumed for the 
prediction of exhauster performance uses flow functions 
with polytropic exponents to describe the flow in four 
regions of the exhauster; namely, nozzle entrance to nozzle 
theoat (n,), nozzle entrance to nozzle exit (nj). secondary 
air from the hood to the annulus at the start of the mixing 
pipe (ng), and the compression process in the diffuser 


(n The relationship between polytropic efficiencies 


ou 


( Uh ) and polytropic exponents (n),for an expansion pro- 


cess,is 
= o 2 Y lous 
Y, = = oe 
~ 8 (x1) %e es Ne 


and for a compression process 


“he —— %. = (ah, =) 











42 


2 
Eckert 2 assumed polytropic efficiencies to be: 


1), = 0.975, ¥) = 0.925, 1) = 0.90, and. ¥) = 0.70. 
Cc 


Secondary Flow 
Tests were conducted to determine ng and Y) q for 
the secondary flow. Secondary flow rate (w,) was maintained 
at 4 lbm/sec, and drive nozzle total pressure ratio 
PIN/P... Was set at 2365) 2745, 2mile;9 1-90, "and ice, Gare 


flow function for this expansion is 


+/ 
wx Vie | 8 = b - BE] (Ee VM _ (6 arg (19) 
Air ire K D le Fro / Fp 


Quantities on the left were either measured or are con- 
stants. The pressures on the right were measured. Ng was 
determined by iteration in computer program PDLY. Results 
of the tests are given in Table I. For 4 lbm/sec secondary 
flow rate,the average values are hg = 1.32 and th = 0.85 
for the supersonic operating range of the exhauster. Con- 
ditions during the test were stable. 

Initial testing to determine ng for secondary flow 
rate of 2 lbm/sec was not immediately successful. In two 


separate tests, results showed polytropic exponents from 


1.08 up to 1.41 and one value of 1.61 which is impossible 


2Onckert, op. cit., p. 98% 
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since ng cannot be larger than &Y = 1.41. During the 
test, annulus static pressure was unsteady and an average 
reading was used for data reduction. Figure 13 indicates 
that an unusual condition exists in the mixing pipe between 
pressure ratios of 2.45 and 2.18. The curves showing mix- 
ing pipe static pressures cross each other for the given 
conditions. This phenomenon repeated itself on four sepa- 
rate occasions. 

In order to determine the effect of a small error in 
the annulus static pressure readings, the computer program 
POLY was used. Values for hood total temperature, total 
pressure, and secondary flow rate were held constant and 
values of annulus static pressure varied in increments of 
O.1 inch mercury or 0.05 psia. The result is given in 
Table II. It is evident that annulus static pressure is 
critical to the calculation of ng since a change of 
*. Oe 1 psia from the test value changes ng between 1.39 
and 1.24. In the final test it was attempted to take all 
readings simultaneously. Results of this test are con- 
sidered accurate and are given in Table III. They show 
that ng, is about 1.28, with 1 4g equal to 0.84. The 
tests showed that the polytropic efficiency was independent 
of drive nozzle total pressure but dependent on secondary 


flow rate. 


Diffuser 
Diffuser polytropic efficiency was determined empiri- 


cally from the static pressure rise between the mixing pipe 


AA 


station M9 to the diffuser exit for all operating conditions 


tested. Figure 27 is a graph of the diffuser static pres- 


sure rise for secondary flow rates of 4 and 2 lbm/sec for 


different pressure ratios of drive nozzle total pressure and 
atmospheric pressure. 

Initially, the exhauster performance prediction program 
was used with the actual exhauster operating conditions as 
inputs and with a diffuser efficiency Nn. = 0.7 as assumed 
by Eckert. The predicted static pressure rise from M7 to 
atmospheric conditions based on Y), of 0.70 is plotted also 
in Fig. 27 for 4 and 2 lbm/sec secondary flow, showing that 
the diffuser efficiency is larger than 0.7. A diffuser 
efficiency of 0.90 was assumed, and the results of the cal- 
culated static pressure rise in the diffuser is shown in 
Fig. 28. Agreement between calculated and actual diffuser 
static pressure rise is good, but the graph shows that there 
is a small variation of diffuser efficiency with drive noz- 
Zzle total pressure ratio. A diffuser efficiency of 0.90 
very accurately predicts the performance at the lower drive 
nozzle total pressure ratios up to about 2.1 but under- 
estimates the performance at the higher pressure ratios. 

The prediction program was altered to vary //, with drive 
nozZle total pressure. Values used for diffuser efficiency 
at different drive nozzle total pressure ratios are shown 
in Table IV. The results of calculated pressure rise in 
the diffuser compared to test data are in Fig. 28. The 


curves show good agreement between calculated and test 


results. 
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It was noted also that the diffuser efficiency varied 
slightly with secondary flow rate. This effect is small 
and need not be taken into account in the simplified model 


used for the prediction program. 


9. Determination of Average Shear and Friction Coefficient. 
Mathematical solutions do not exist that exactly pre- 
dict the conditions ina flow as complicated as that in the 
mixing tube of the exhauster. However, it was attempted to 
establish a value for the average shear stresses in the 
flow by an approximate method. Such an approximation was 
obtained by not attempting to satisfy the hydrodynamic dif- 
ferential equations for each particle but by integrating 
the profile distributions obtained in laboratory tests and 
hence solving the momentum equation for the flow between 


two stations. The momentum equation 


eM, + [+ dA, = |v, dyn, + | aa, <A, (20) 
£ ® 


0) @ 
where: 1 represents conditions at the upstream station 
2 represents conditions at the downstream station 


eS surface area of the mixing pipe between 
stations 


can be used, but in turbulent flow with large velocity gra- 
dients oe consists of the "true shear" due to the visco- 
sity of the fluid and the "apparent shear" due to the losses 


produced by the turbulent interchange of fluid particles. 
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Prandtl states~/ that except in layers directly next to the 


wall, the true shear is small enough to be neglected in com- 
parison to the apparent shear. ‘The approach assumes axi- 
symmetric flow, a perfect gas, axial velocities everywhere 
and steady flow, with gravitational effects being ignored. 
Equation 20 is used for the profiles of velocity, pressure 
(total and static), and temperature measured at stations 

Mi, M4, and M8. 


With the pressures in psi and r in inches, there are 


ABA = 2arrdr 


M4 yr dy . aw / PV rdr 
din = PVIA = GP Vata te « 22 (> + ) 


Vdm = a (rar tvirdy | 


Equation 20 becomes 








2 2 pA 
oa une rar + ik rdr = a | ad ne” el el rdy + T A. 
. 2. 
O @ 


Numerical integration of large scale reproductions of Figs. 


29 and 30 gave at station Ml 


+ g | AM a + | trae = 342.5 +475.0 = 817.5 Ibs. 
© O 


2 le ded L., The Mechanics of Viscous Fluids, 
Aerodynamic Theory (California: Durand Reprinting 


Committee, California Institute of Technology, January 
745), VOl. III G, p. 135. 


‘AT 





2 
ATT ) AO dy i f rd = 3344 + S525,5 = 869.9 bs, 
@ " @ 


and at M& 


2 
a fs Ve dy + R ray = 222.0 + S465.3= 787,43 lee 
Tg @ 


The results, which show an increase in momentum from 


station Ml to M4, must be in error although use of the 
momentum equation to model the mixing pipe performance and 
the mathematical integration can be considered correct. 

The error could have been introduced by the blockage effect 
of the traverse probe and the variation of the overall 
exhauster performance with probe position. To illustrate 
this condition, Table V shows the values of the hood total 
pressure for the same exhauster operating conditions, but 
with the traverse probe at different locations. A varia- 


tion of 1.65 in. Hg. hood total pressure occurs between 


the conditions where the traverse probe is at station Ml, 
and the probe is completely removed. Another possible 
reason for the error is that it was not possible to produce 
exactly the same settings of pressures and temperatures for 
each traverse since, to change traverse probe locations, it 
was necessary to shut down the turbine and exhauster. Re- 
establishment of all previous conditions was extremely 


difficult. Analysis of the data indicates, however, that 


this effect should be secondary. 
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The calculation was continued between station M4 and 


M8 using the integrated values obtained. 


F = 869.9- 1787.3 = 82.6 Ibs. 


RA. = 6.46 cota 


m4, - Me 


T = 8.19 8/4? 


From the test data and results obtained, the turbulent 
friction coefficient f£ was also calculated. Using 

28 
Vavra, the pressure drop QO p due to friction may be 


approximated by 














Bae 
Bie 
AP = + pe (psf) 
where f = turbulent friction coefficient (dimension- 
less) 
L = length of pipe (in.) 
D = diameter of pipe (in.) 
P = average density between stations (slugs/£t?) 
V = average velocity between stations (ft/sec) 
with 
* (HY) _  ¢.38 ® ( >) 
P= RT 9 = 60 + Slugs /$T 
": a1 oat 
= -t 2 
Va 


SU auica, M. H., Theoretical Evaluation of Exhauster 
System of Transonic Turbine Test Rig (TNO5T, U. S. Naval 
Postgraduate School, Monterey, California, 1965), p. 26. 
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The total frictional force may be expressed in terms of a 
pressure drop 


& 
FO = i AP =0.303 ap (\bs.) 


With I determined, £ becomes 


T_F, 


|p Ve (25 


£ = 585 
Average values used to calculate f£ were: Gy 530°R, 


p = 12.4 psia, and V = 700 ft/sec. 
f = 0.106 


This is one order of magnitude larger than the friction 
factor 0.017 predicted by the Moody diagrance for commer- 
Cial steel pipes. The value should have been closer even 
with velocities of 700 ft/sec. Shapiro cites numerous 


30 that confirm the use of the incompressible friction 


tests 
formula for conditions up to Mach = 1.0. 
The calculations were repeated between Ml and M8. 
With F, = 29.7 lbs., L = 86 in., D = 11.3 psia, T = 500°R, 
and V = 750 ft/sec, Eq. 21 gives 
f = 0.024 
This value is still large but has the correct order of 


magnitude for the friction coefficient. 


-2shapiro, op. cit., Fig. 6.15, p. 184. 





30-n34., p. 1071. 
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The results of the calculations for friction coeffi- 
cient suggest that the traverse data at station M4 are not 
accurate. Better results were obtained between Ml and M8 
than M4 and M8. No particular quantity in the M4 data was 
considered unreasonable. The conclusion is that the cal- 
culation of an average shear force and friction coefficient 
for the mixing pipe,using the ejector test data, produced 


unreliable results. 


10. Evaluation and Discussion of Predicted Results. 


Actual test conditions were input to the exhauster 
prediction program to determine a predicted hood total 
pressure. Three predictions were made: Program EJECT 1 
used constant polytropic efficiencies determined from test 
data for the secondary air and diffuser and an assumed 
drive nozzle polytropic efficiency of 0.975, EJECT 2 used 
Eckert's polytropic efficiencies, and EJECT 3 used EJECT 1 
efficiencies, except that the diffuser efficiency was 
varied with drive nozzle total pressure from 0.90 to 0.93. 
The results of the predictions for drive nozzle I anda 
comparison with the test results of run 18 are shown in 
Fig. 31 for 4 lbm/sec secondary flow, and in Fig. 32 for 
2 lbm/sec. Appropriate changes were made in the program 
to reflect the new dimensions of drive nozzle II. The 
prediction compared to test results is shown.in Fig. 32. 

The results show that Eckert's assumed polytropic 


efficiencies greatly underestimate the exhauster performance. 
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At high drive nozzle pressure ratios, the use of the theo- 
‘retical prediction program EJECT I with experimentally 
determined exponents very closely establishes the actual 
performance of the exhauster. The accuracy of the predic- 
tion decreases as drive nozzle pressure ratio decreases. 
At the low pressure ratios the exhauster performs much 
better than predicted. 

Much of the difference noted above can be attributed 
to the physical formation of an oblique shock within the 
drive nozzle for off-design operation. The classical 
approach to describing nozzle operation notes a third 
critical pressure which is the exit static pressure that 
provides design operation. As the receiver static pres- 
sure is increased, an oblique shock forms, originating 
at the exit plane edge. Further increase in static pressure 
increases the angle of the oblique shock until a normal 
shock forms in the exit plane for the second critical pres- 
sure. Still further increase in pressure causes the normal 


shock to move from the exit toward the nozzle throat. This 
idealized flow development does not occur in actuality 


because of friction and three-dimensional flow effects. 

The latter, in particular, produce: internal oblique shocks 
which have the effect of extending the shock action Over an 
extended axial distance in the nozzle. Tests conducted by 
Stodola>* show internal shock effects over the axial distance 


of a nozzle. 


Sliall, N. A., Thermodynamics of Fluid Flow (Englewood 
Chip Sec e Niwas Prontice Hall. ine., 1951), Fig. 9.5, paige 
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By removing the mixing pipe and operating as a free 
jet, the presence of an internal oblique shock in the drive 
nozZle could be made visible because of condensations of 
the atmospheric humidity in the primary stream. Figure 33 
is a photo of the free jet operation. The oblique shocks 
originated inside the nozzle and extended outside the 
nozzle to their intersection with a near normal circular 
shock in the center area of the nozzle. A sketch of the 
observed pattern is shown in Fig. 34. With increased drive 
nozzle pressure, the entire shock structure moved towards 
the nozzle exit. 

To determine the displacement of this shock structure 
in the drive nozzle under exhauster test conditions, six 
static pressure taps were arranged in drive nozzle I in 
positions shown by Dl through D6 in Fig. 35. At a secon- 
dary flow of 4 lbm/sec the drive nozzle total pressure was 
set at values of PTN/P,,, Of 2.65, 2.45, 2.18, 1.90, and 
1.63 to duplicate conditions of run 18. The measured 
static pressures at the nozzle walls are shown in Fig. 36. 
They show that the internal shock moves toward the nozzle 
exit if the pressure ratio is increased. It should be 
noted that for each condition the pressure at the discharge 
of the secondary nozzle was below that required for oblique 
shock formation at the exit plane for isentropic flow 
theory; that is, the pressure was always less than the gnd 


critical pressure. 


ao) 


During the free jet operation it was observed also that 
the primary jet did not completely fill the exit plane when 
Oblique shocks occurred inside the nozzle. There existed a 
slowed down region of flow at the outer perimeter of the 
nozzle discharge that had a radial thickness of 0.5 to 0.75 
inch. This phenomenon can be noticed in Fig. 34 and must 
be responsible for the apparent thickness of the mixing zone 
at the nozzle exit which is illustrated in Fig. 24. 

The order of magnitude of expected boundary layer ~ 
thickness at the nozzle exit was determined with Eq. 22 
which can be used to evaluate the boundary layer thickness 


32 


on a flat plate in turbulent flow. The Reynolds number 


applied is based on distance from the nozzle throat to its 


exit. 
$2 oan(utys - 
h oe 
$§ = boundary layer thickness 
iw = kinematic viscosity (ft/sec?) 
= distance from throat to exit (ft) 
V = velocity (ft/sec) 


Based on average conditions 


§ © 0.07 inch 
This boundary layer thickness igs much smaller than the ob- 


served distance of about 0.5 to 0.75 inch. The increased 





32gcnlichting, H., Boundary Layer Theory (New York: 
McGraw Hill Book Co. Inc., 1955), Eq. 2.9, p. 38. 
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boundary layer thickness must be caused by flow separations 
from the nozzle walls. This separation is the result of 
the internal oblique shock structure that exists in the 
drive nozzle. The nearer the throat the separation occurs, 
the less of the diverging portion of the nozzle is used to 
increase Mach number of the flow; that is, full expansion 
is not realized. The effect is readily apparent by ob- 
serving average nozzle Mach numbers obtained at station Ml 
for Pa ee. of 2.65, 2.18, and 1.63. Test data showed 
Mach numbers of 2.08, 1.18, and0O.88 respectively. Figure 37 
is a sketch of the shock conditions that appear to exist. 

The effect of the internal oblique shock and resultant 
primary flow separation is to raise the static pressure of 
a portion of the nozzle exit area to the pressure of the 
static pressure at the discharge of the secondary nozzle. 
The momentum theorem expressed by Eq. 17 uses in the 


actual prediction program 


PSUM Ay = Pin 1 + Pin Aig 
a = mixing pipe area 
Ain = secondary flow annulus area 
A, = nozzle exit area 


Pin = nozzle exit plane static pressure 


Pim = secondary air static pressure. 


In the program the value of PSUM is increased by itera- 
tion until the conditions at the diffuser exit match the 
atmospheric pressure. In the prediction program the 


assumption is made that Py and Pjm are constant, 


N 


She 


respectively, over A, and Aj, in reality P,, occurs 
On an area smaller than Ay. and Pim On an area greater 
than Aim for conditions where the internal oblique shock 
exists. The result is that for a given value of PSUM Ao, 
the pressure Pin may be smaller than predicted to satisfy 
the equality. Thus, the hood total pressure would also be 
lower. This condition should be taken care of in the pro- 


gram for a more correct performance evaluation at low 


pressure ratios. 


11. Conclusions and Recommendations. 

The performance analysis of the exhauster system based 
on polytropic efficiencies and conservation of momentum 
very closely predicts the exhauster performance with re- 
gards to hood total pressure at high drive nozzle pressure 
ratios. The prediction under-estimates exhauster perform- 
ance at the low pressure ratios because of the presence of 
internal oblique shocks in the drive nozzle. Accuracy ob- 
tained by the prediction program was the same for drive 
nozZle II which had dimensions different from nozzle I. 

The hood total pressure is significantly effected by 
the presence of a traverse probe in the mixing pipe. The 
amount of hood pressure rise ig a function of probe posi- 
tion. Therefore, traverse information presented cannot be 
considered as having been obtained for precisely the same 


operating conditions. 
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Traverse profiles near the drive nozzle exhibit steep 
gradients of pressure, temperature, and velocity. Velocity 
profiles obtained appeared irregular. ‘The lower centerline 
velocity is probably due to the shock structure that occurs 
in the nozzle. It appears that a sufficient condition to 
determine the start of the main mixing region in a restric- 
ted jet is to determine where the velocity profile may be 
non-dimensionalized to a single profile. Results from flow 
traverses showed that the initial mixing zone expands 
radially at a faster rate than predicted by Victorin's 
incompressible approximation. 

Secondary flow polytropic efficiencies are independent 
Of drive nozZle pressure but dependent on secondary flow 
rate. Diffuser polytropic efficiency is independent of 
flow rate and slightly dependent on drive nozzle pressure. 
The diffuser polytropic efficiency is about 0.90 which is 
considerably higher than the previous estimate of 0.70. 


Little test data are available on the behavior of com- 
pressible free jets. The facility used for the ejector 
tests could be adapted to obtain test data for experiments 
in this area. In addition, further insight into the shock 
pattern in the drive nozzle and the mixing pipe might be 
gained if an adaptation for Schlieren photography were de- 
Signed. A- study of the effects of the frictional and 
three-~dimensional effects on shock formation in a circular 
nozzle is considered the next step necessary to improve the 
exhauster prediction analysis so that it produces more 


accurate results for the entire operating range. 
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ANNULUS POLYTROPIC EXPONENT 


STATIC PRESSURE Secondary Flow 
(PSIA) (ng) 
6.47 1.47 
6.42 Less 
6.37 1.34 
6.32** 1.30 
6.27 1Lii2a 
6.22 1.24 
Saaky La2z2 


**actual Test Reading 


VARIATION OF Ng 
WITH ASSUMED ANNULUS STATIC PRESSURE 


Wry = 1.92 lbm/sec TTD = 547.5° PTD = 6.8 PSIA 


TABLE II. 
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PTW Poe if... 


2.65 O220 
2.45 Omg: 
2. WS 0.92 
bo20 0.98 
1.63 0.94 


Variation of Diffuser Polytropic Efficiency 


With Drive Nozzle Pressure Ratio 


TABLE IV. 
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Hood Total Traverse 
Pressure Probe Location 


Pas 7 Py ae 


=529 in, Hq. M1 
-7.0 M2 
-7.2 M3 
— Th me. M4 
-7.85 M5 
-8.15 M6 
“8.55 No Probe 


Hood Total Pressure 
versus 


Traverse Probe Location 
PTN = 32 psia We = 4 lbm/sec 


EJECTOR TESTS 


TABLE V. 
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Diffuser 


ixing Pipe 
Turbine Iniet Exhauster 7 Rotor Support 
Manual Valve = Bearing 
tence Shaft Seal 
Dynamometer 
Dynamometer 


Support Bearing 
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APPENDIX A 


Exhauster Performance Analysis 


An analysis was made of the exhauster performance to 
predict hood total pressure which essentially parallels the 
preliminary design calculations of Vavra. However, the 
previous work assumed that the nozzle operates at design 
conditions and that the static pressure at the entrance to 
the mixing section is constant with radius. The analysis 
explained in the following permits a discontinuity in sta- 
tic pressure between the nozzle exit (Area Al) and the 
secondary flow annulus (Area A1T) at the start of the mix- 
ing pipe. | 

In order to simplify the equations, polytropic pro- 
cess relations as well as non-dimensional flow functions 
were used. These were taken from Vavra. 


For an expansion process 


ace 
ane [E - 7 = |e | ae 2)" _(6 \m — 
2 Pry ps Fy 


where 





yx 
_ (A2) 
Ys: -T 
z sly) 
¥- 7, ( 
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n = polytropic exponent 





UE = polytropic efficiency | aT 


aT ss 
% = ratio of specific heats 
w = flow rate (lbm/sec) 
TE, = inlet stagnation temperature (°R) 
13) ee Flow area at station of interest (in.*) 
Py = Inlet stagnation pressure (psia) 
O 
P59 = static pressure at station of interest 
(psia) 
g = universal gravitational constant 8 
lbs sec? 


Subscript 1 refers to inlet conditions 


Subscript 2 refers to station of interest 


For a compression process 








wT \R -$ = 
RA, Jy te i 
where 

P, = static pressure at l 


Y) __« (A4) 
|-2x(i-NQ,) 


The flow through the nozzle will be choked when it is 


operating at supersonic conditions. Therefore So (on 


2 
$. ) must be a maximum value. Differentiating $. and 
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setting the result equal to zero yield the critical pres- 


sure ratio for choked conditions 


v1 


R = (2 \"@ 
lm N+i 


eeu 





Substituting this value into Eq. Al,the expression for the 
flow function of a choked nozzle simplifies to 


n-! 
$ = Fa 2% (2 ; (AS) 

choked N+l x) Via 

Nozzle 

It was assumed that the drive nozzle total temperature 
and total pressure were known,as well as the secondary flow 
rate and total temperature. Polytropic efficiencies for 
the drive nozzle are assumed or determined from testing and 
held constant over the operating range of the exhauster. 


Notation in the development is consistent with Fig. 3. 


Drive Nozzle 


For a choked nozzle, 
| 


$= (4 a 2x (ais t (a6) 
N N+ | x-{ tee 


The flow functions from the nozzle entrance to the throat 





( Q. ), and from the nozzle entrance to the exit ( oes) are, 


respectively, 
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¢ _ Wy View ing 
a) 


A, Ry 3 ale 
$= Ya Vw ye 
Cx A, 3. 9 (A8) 


There are 


‘a A. 5 eA, (A9) 
and 
Pe, = $= (A10) 


Expressing the exit flow function in terms of the pressure 


ratio Poy’? rN 
£..* Ver (8): GE) 
ex oll a \ale om 


With Foy given by Eq. AlO, and with PTN known, PIN can be 





solved by iteration. Then 


tT Fi oe (A12) 

IN Tony Peng 

VIN 7 V2aJde, (Tn-Tin) i 
; — Va i ) (A1l4) 
Wy 7 ie En A, oe 
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Secondary Flow 


Secondary air is that part of the high pressure air 
provided by the Allis-~Chalmers compressor that passes 
through the test turbine to the hood. For the expansion 


process from the hood to the annulus area 
ar we “Ing Pr Poy [/Pe Vm, /Poltatld - 
a Ne (A15) 


With the secondary flow rate and hood total temperatures 
Known, a value of PIT was assumed, and PTD is calculated 


by iteration and matching. For the assumed PIT 


ny ol 
Ma 








(A16) 


2 Kir 
1% . >. ; 
TO 


V. = V2qJc, bs Gemelices (A17) 


Mixing Pipe 


The momentum equation was used for the mixing pipe to 
obtain conditions at the diffuser entrance. The flow rate 
through the drive nozzle is known from preceding calcula- 
tions and values must be assumed for the secondary flow 
conditions. Applying the momentum equation to the control 
volume between the nozzle exit and annulus plane (sta- 


tion 1), and the entrance to the diffuser (station 2), 
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— Fe (A18) 


r 


Wry +W Ww | 
(atuttaly, an — aN, = Fr Aly My a BA 


Equation Al&8 holds for uniform velocity profiles and uniform 
static pressure distribution over each area. 
The friction force Fr may be expressed as a pressure 


loss 
TT 2 
Fe = — D, Ap = af fl. (A19) 


where 


(A20) 


rol 


AP = f= 


IL = length of duct 
D5 = duct diameter 

V = average velocity 
P = average density 
f = friction factor 


The friction factor was determined to be 0.017, in the man- 
ner outlined by Naviaux. >? The friction force was figured 
with one iteration only. For the initial solution of 
Eq. Al8, Fe was considered zero. After obtaining the con- 
ditions at the diffuser entrance, average values of Vand 
pP for the mixing pipe were calculated,and a AP found 
using Eq. A20. AP Az was then subtracted from P,,)Aj,, + 
PinAjy to obtain a new solution of Eq. Al18g. 

33Naviaux, op. cit., pp. 25 =— 26. 
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Unknown quantities in Eq. Al18 are Vg and py. The 


pressure p5 can be expressed by Vy to obtain one inde- 


pendent variable if the total temperature after mixing is 


taken as the mass weighted average of 


a W, Trp te Wey eas 
Te 7 (A21) 


Wy + We 


This condition assumes the existence of an adiabatic pro- 


cess. With 








we + Ww, = W (A22) 
and the continuity equation 
W = P A. V, (A23) 
, : V a 
y= wR Ip _ wR Tr ——_ 
2. A, Bb Ae PA 241 Cp (A24) 
Therefore 
: 2 
fs = wh | Tre — af (A25) 
Substituting Eq. A25 into Eq. Als 
were wy My =p a tea - WRI ve (426) 
g 4 IN gq fit IT CIT NTT V2 T2 24T cp 


LO® 


Multiply Eq. A26 by V5 and rearrange 


F-a8 | ; fir, + Be, er aA), *WRip2 = O (A27) 








Let 
ee 
ey ' 
Bo - -|*V, met RPA PA (A29) 
p= WRE, (A30) 
By v5 + By V, + B3 = 0 (A31) 
Thus 
ee Se (#2 Bs (asa) 
Pa) RB, 28, B, 
and 
v4 
Es — Vai 
Th deg 294Jcp = 
a 
2 dike | 
Ae 7 a (7 — 


The negative sign before the square root yields the 
subsonic solution. Supersonic conditions at the diffuser 


are outside to operating range of the TTTR installation. 
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Finally, for the compression in the subsonic diffuser 
RON yo wil ower ee 
-te c=) Sa RM Pr, n+ 
P ve (A35) 
A, ry 1 al 


P. is solved by iteration and compared to atmospheric 








pressure. If it does not match for the chosen values of 
secondary flow conditions, a new PIT is selected and the 
process repeated until P3 equals Potm: With pz 


obtained 


n.-! 


f \™ 
oi z Te (2 | : (A36) 


< 


aqJep (Ty, —Thy} _ 


ie es 
i = f (—* 3 
T3 .\ 7, (38) 
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APPENDIX B 
Exhauster Performance Prediction Program 


The analysis of exhauster performance presented in 
detail in Appendix A was programed in Fortran IV for cal- 
culation on the IBM 360 computer. Input data include the 
assumed polytropic efficiencies, drive nozzle total pres- 
sure, drive nozzle total temperature, secondary air total 
temperature, secondary flow rate and atmospheric pressure. 
The principal quantity to be determined is the predicted 
hood total pressure. The printed output does include 
final values of temperature, velocity and pressure calcu- 
lated during the computations. 

The program uses an executive MAIN routine to arrange 
the sequence of calculations, but all calculations are 
performed in subroutines. The name and function of each 
subroutine is: 

CSTANT; a grouping of all the basic quantities 
used in the entire program. 

CHOPHI; a calculation of the non-dimensional flow 
function for a choked nozzle. 

NOZZLE; computes conditions at the nozzle exit. 

SECFLO; computes conditions at the secondary flow 
annulus area at the entrance to the mixing pipe. 

MMENTM; Using conditions determined by NOZZLE and 


SECFLO, conditions at the diffuser entrance are calculated 


Tike 


based on conservation of momentum. One iteration is in- 
cluded to determine the wall frictional shear force, 
SBCOMP; computes the compression process in the 


diffuser to obtain conditions at the diffuser exit. 


Input information required to operate the prediction 
program is: 
Gard 1; 


ETAD - Secondary flow polytropic efficiency 


ETAN - Nozzle polytropic efficiency, entrance to 
throat 
ETAL - Nozzle polytropic efficiency, entrance to 
exit 
NSETS ~ Number of sets of data to be computed 
Garda 2; 


PTN - Drive nozzle total pressure (psia) 


TIN -— Drive nozzle total temperature (°R) 

TTD - Hood (secondary air) total temperature (°R) 
WT - Secondary flow rate 

PA ~- Atmospheric pressure 


Program EJECTOR with a sample output sheet is given on the 


following pages. 
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APPENDIX C 
Data Reduction Programs 


Data obtained from laboratory tests to compute secon- 
dary flow rate, velocity profiles, and polytropic exponents 
were reduced by the computer programs included in this 
section. Input requirements are listed on comment cards at 
the start of each program. 

Program SECFLO computes the mass of air that passes 
from the hood to the mixing pipe. The program basically 
calculates Eq. 1 of the main text. In addition, the small 
quantity of air that enters the hood through the shaft 
labyrinth is approximated by Eq. 2. Naviaux's flow nozzle 
coefficient is determined as a function of Reynolds number 
by a series of "IF" statements. The values used range from 
1.00 - 1.05 in increments of .Ol. 

Program VELOCITY computes Mach number and Velocity 
from traverse probe data. Input pressures are tested for 
Supersonic or subsonic conditions. Either Eq. 4 or Eq. 5 
is used depending on the condition that exists. 

Program POLY computes by iteration the secondary flow 
polytropic exponent from test data. Eq. 19 with all quan- 


tities known but nN, was used in the calculation. 
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APPENDIX D 
Restricted Jet Operation 


In order that the reader might better understand the 
general terms used in this thesis to describe jet operation 
and associated mixing phenomena, the following general 
description of jet operation is presented. 

Figure Dl of this appendix pictures the physical sit- 
uation to be discussed. Essentially, two co-flowing 
streams enter a mixing pipe. Abramovich’~ states that with 
this situation, the velocity components along the y-axis 
(the radial components for the exhauster) prove so small 
that they can be disregarded in engineering problems in- 
volving jet theory. 

The two co-flowing streams enter the mixing pipe with 
their fluid particles traveling at different velocities. 
The viscous forces acting between the two flows cause a 
slowing down of the particles of the primary jet and a 
speeding up of the secondary flow particles. The area of 
the flow where the interaction between the particles of 
the two streams takes place is called the mixing zone. 

The mixing zone thickens with axial distance, causing a 
gradual “eating up" of the primary and secondary core 
regions. 

Two main mixing regions are designated in the mixing 


pipe. The initial region hasthe axial distance where the 


33apramovich, @P.wei tsp.» 5D 
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secondary core region and the primary core region still 
exist. The main mixing region is that part of the pipe 
where the mixing zone has touched the pipe walls and the 
primary core region no longer exists. The regions are 
shown in Fig. Dl. In the literature use is sometimes made 
of a transition zone which is the region that occurs if the 
mixing zone contacts the pipe wall before the primary core 
is "eaten up" or vice versa. This thesis uses the simpli- 
fied model in which it is assumed that the length of the 
transitional region is equal to zero. In the main mixing 
region, the jet becomes similar in appearance to a flow 
of fluid from a source of infinitely small thickness (in 


an axially symmetric case,the source is a point). 
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pressure ratios available for turbine testing. Exhauster 
operation results in interaction between two co-flowing streams 
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subsonic. This thesis is an experimental study of the mixing 
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were compared to test results. 
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